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Peo RACT 


An experimental investigation was conducted to examine 
convection heat transfer in a duct of rectangular cross-sec- 
tion, having both a straight and curved flow passage. One 
wall was kept at a constant temperature, with the opposite 
wall being adiabatic. Heating of the air flowing in the 
channel was accomplished through Joulean heating of a wall 
composed of Temsheet. 

The experiments were conducted for steady state in 
Laminar and Transition flows. Development and progagation of 
the Taylor-Gortler vortices was shown to enhance the heat 
transfer rate in the curved section as compared to the 


straight section. 
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I. INTRODUCTION 


meme y LOR-GORTLER VORTICES 

cmce the carly part Of this century, considerable data 
and numerous studies have demonstrated that fully developed 
laminar flow along a concave wall does not remain two-dimen- 
fwenmal |[Rets. 1, Z, 3j. Instead, the flow forms a system of 
memeealeyOrtices, of Counter rotating pairs, whose axes are 
Paeeom@ed in the direction of the principle fluid flow. This 
phenomenon, Known as Taylor-Gortler vortices, is the result 
of the variations in the centifugal forces acting on the 
fluid particles. Schematics of this type of motion are shown 
meereures |] and Zz. 

In a channel that is curved in the streamwise direction, 
the fluid particles near the center of the flow cross-section 
are subjected to higher centifugal forces than those slower 
moving particles nearer to the boundary wall. Consequently, 
the tendency is for the fluid in or near the center of the 
Channel to move outwardly towards the concave wall. This 
action then causes the fluid particles near the boundary wall 
to move in the spanwise direction, and finally move radially 
inward, replacing the outwardly moving particles. Once in the 
center of the channel these fluid particles come under those 
same stronger forces of the free stream velocity and continue 


the rotational effect. This cyclic motion causes the 
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Figure 2. Schematic of Taylor-Gortler vortices between 
parallel plates 








[ieitetenweanad propagation, of the counter-rotating Taylor- 
Gortler vortices. 

PeenouCMmeiemlaylor-Gortler vortices are primarily a 
laminar flow phenomenon they may affect the transition from 
laminar to turbulent flows, and are therefore important to 
understand [Ref. 2]. The additional mixing provided by the 
secondary motion has been thought to account for the increase 
mn the rate of heat transfer from concave curved walls versus 
Someht walls [Ref. 4]. 

The Taylor-Gortler vortices appear to have similarities 
memomener vOrtex flow patterns such as the longitudinal vor- 
tex rolls developed in the laminar forced convection heating 
of fluid layers between parallel plates [Rem |, it has 
also been thought that the cross-hatching, observed in 
reentry vehicles, can be explained at least in part, by the 
presence of streamwise vortices [Ref. 6]. 

There are many possible applications that could result 
from a more thorough understanding of the Taylor-Gortler 
vortices and their effect on heat transfer and fluid flow 
characteristics. Applications, such as improved turbine 
Oieee Cooling {Refs. 7, 8, 9], heat exchanger designs, as well 
as other such engineering designs and capabilities could take 
advantage of these vortices and improved heat transfer 


cChawdeteristresmof the flow in the presence of these vortices. 


I 








Breads LORY. 

The instability of an inviscid fluid rotating symmetri- 
Smebyeabout an axis, aS in Curved flow, was first considered 
Pyeeeerd Rayleigh in 1916 eet. 10} In assuming that the 
mmo was NnOn-viscous, he determined that a necessary 
Sememcton, tor the stability of an inviscid fluid motion, was 
maeemene Circulation increased with increasing radius. G.I. 
ieee r, in 1923, [Refs. 1, ll] expanded the earlier works of 
bey eetmeoh by an extensive analytical and experimental study of 
viscous fluids. His investigations focused on the flow 
between two cylinders, in which the inner cylinder rotated 
while the outer cylinder remained stationary. Taylor 
ascertained that such couette flows became unstable when the 


value of the dimensionless Taylor number, defined as: 


Ta 





exceeded a critical value of 41.3. In the above equation, 
'd' is the width of the gap, assumed small when compared to 
‘Ry ', Pacerad@ius of the inner cylinder, ce is the periphal 
velocity of the inner cylinder, and 'Re' is the Reynolds 
number. Taylor determined that for those cases in which the 
Vaalewoeethe Taylor number became greater than the critical 


value, a secondary motion developes and the Taylor vortices 
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Poe ittyvewOrea similiar nature is also observed when a 
viscous fluid flows in a curved channel due to a pressure 
gradient acting along the channel wall. This problem was 
considered analytically first by W. R. Dean [Ref. 12] in 
mozG, for a channel formed by two concentric cylinders, 
where the radius of the inner cylinder was large in comparison 
to the small spacing between the inner and outer cylinder 
walls. Dean, in his studies, concluded that there would be 
an initiation of the flow instability, and the formation of 
vortices, (similiar to the Taylor vortices), when the Dean 


number, defined as: 


d 
De = Re 


R. 
iL 


memmevead 2 Value greater than 36. In the above equation, the 
Reynolds number is based on the mean velocity of the unper- 
P@eoed tlow, ‘d" represents the channel half-width, and 'R;' 
1s again the inner cylinder radius. The analytical work of 
Dean, was later verified by W. H. Reid [Ref. 13], using an 
approximate numerical solution. 

H. Gortler [Ref. 2], in 1940, studied the stability of 
laminar boundary layer profiles, on curved walls under the 
Meemence Of small disturbances. In his study, Gortler found 
that these disturbances were similiar in nature to those 
studied by G. I. Tavlor, which had led to the fluid flow 


instability in the form of vortices. By approximate numerical 


ley 
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Ss wremlattons, Gortler concluded that only on the concave 
boundary walls were the amplified disturbances produced, and 
that the overall flow profile appeared to remain laminar in 
Mature. In 1955, Gortler's approximate results were verified, 
teen exact solution, by G. Hammerlin, as reported by H. 
Schlichting [Ref. 14] and then further substantiated in a 
more extensive numerical analysis by A. M. O. Smith [Ref. 3]. 
More recently, verification of these numerical solutions have 
been demonstrated by use of hot wire anemometry, laser doppler 
Systems, and flow visualization techniques, [Refs. 15, 16], 
and in 1976, Y. Aihara [Ref. 17] conducted a non-linear 
analysis of the Gortler vortices. 

With the growing interest in the effects of the secondary 
flows associated with the Taylor-Gortler vortices, numerous 
specific studies have been published concerning the influences 
of these vortices on the transfer of heat in both the laminar 
and turbulent flow regimes. F. Kreith [Ref. 4], in 1955, 
studied the influence of heat transfer with respect to the 
curvature of the boundary wall for fully turbulent flows, 
and concluded that the heat transfer from the heated concave 
boundary wall was considerably higher than that transferred 
from the convex boundary wall of the same curvature and under 
Pim@liar turbulent flow conditions. 

medion, Lf. rersen [Ref. 18], considering the special 
cases of very high and very low Prandtl number fluids, 


related the increase in heat transfer rate from a curved wall 
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to the presence of the Taylor-Gortler vortices. His con- 
Glusions were that the effect, of the development of the 
Mmemerees, Was tO increase the rate of heat transfer through 
the boundary layer. 

One of the first experimental works on the effects of the 
Taylor-Gortler vortices with respect to enhanced heat transfer 
was by P. McCormack, et al., [Ref. 19], published in 1979. 
With the data, they concluded that it was necessary to retain 
the non-linear terms of the flow equation in order that the 
complete theoretical solution of the effects of the Taylor- 
Gortler vortices on heat transfer be explained. In 1977, 

R. Kahawita and R. Meroney [Ref. 20], furthered McCormack's 
work and concluded that the higher order terms, as well as 
the normal velocity components of the main flow, became 
miemersingly more Significant to the calculations at smaller 
wave numbers. 

In considering the forced convective heat transfer for 
fully turbulent flows, between horizonal parallel flat plates, 
Mier ae ana Y. Uchida [Ref. 5], in 1966, reported the 
formation of longitudinal vortex rolls with an axis parallel 
to the flow direction. These streamwise-axis directional 
rolls, similiar to the Taylor-Gortler vortices, formed when 
a critical value of temperature difference between the 
plates was exceeded. In 1971, M. Akiyama, et al., [Ref. 21], 


Sonrrmed tne formation of these vortices at a critical 


Ee) 








Rayleigh number of 1708, for the case of heating from the 
bottom plate surface. 

Siemewiacmpeen a rather limited amount of published 
literature regarding the flow and heat transfer in curved 
Smemanels Of rectangular cross-section. Much of the pub- 
lished literature involves the development of numerical 
approximations and solutions for heat and mass transfer in 
curved ducts of various geometries (most notably, that of 
Seeeetmelar and triangular ducts). K. Cheng and M. Akiyama 
[Ref. 22] developed a numerical solution for forced 
convection heat transfer with laminar flows in curved channels 
Sewer ectangular cross-section, but only for small aspect 
ratios, (B/D; where, 'B' is the channel width and 'D' is the 
G@emen of tluid flow). Using the fully elliptical forms of 
the transport equations, G. Yee and J. A. C. Humphrey 
{[Ref. 23], developed a numerical solution for flows and 
boundary conditions similiar to Cheng and Akiyama. kK. Cheng, 
et al., [Ref. 24], in 1974, using numerical methods solved the 
Graetz problem for a curved channel of square cross-section 
(B/D=1). 

meoooe >. SeePetukhoy and V. N. Popov [Ref. 25], obtained 
an analytical expression for the Nusselt number for the flow 
Of an incompressible fluid with variable properties ina 
@eecular tube. In 1976, A. A. Shibani and M. N. Ozisik 


[Ref. 26] using matched asymptotic expansion techniques for 
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the case of uniform wall temperature, solved the heat 
transfer problem between parallel plates with turbulent flow, 
for a wide range of Prandtl numbers. 

As to analytical and experimental studies of the parallel 
fees channel flow heat transfer problem, Y. Mori, et al., 
fen? ), Obtained results for hydrodynamically fully 
developed flows with constant wall heat flux, in curved 
channels of square cross-section. W. M. Kays and E. Y. 
Leung, [Ref. 28], experimentally obtained solutions for 
Mmiemimeent tlow heat transfer in a concentric circular tube 
annulus with fully developed velocity profile and constant 
Memmi dtce per unit length, for a fluid of Prandtl number 0.7. 
feeeeeeeao |Ret. 29] and J. Ballard [Ref. 30], reported results 
for large aspect ratio (B/D = 40) rectangular cross-sectional 
Channels with laminar flows. P. F. Brinich and R. W. Graham 
[Ref. 31], reported results for turbulent flows in a rectang- 
ular curved channel with an aspect ratio of 6, for the inner 
wall heated, the outer wall heated, and both walls heated. 

In 1979, F. D. Haynes and G. D. Ashton [Ref. 32], in studying 
heat transfer from a river to its ice cover, reported results 
for rectangular cross-section channels of aspect ratio 10, 
and compared them to the data obtained by Ashton [Ref. 33], 


and Hsu [Ref. 34]. 
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Pein. OF TERE STUDY 


The primary purpose of this study was to investigate the 
effect of the Taylor-Gortler vortices, on the enhancement of 
@iem@eat transier rate in a curved channel of rectangular 
cross-section and large aspect ratio. This investigation 
began with flow velocities that were in the laminar flow 
regime and continued into the transition region, between 
laminar and turbulent flow. Further, this study compared the 
heat ransfer rate obtained in a curved section with the rate 
@btained in a straight section of identical aspect ratio. 

It was expected that the presence of the Taylor-Gortler 
vortices would enhance the heat transfer process at all 
velocities in the flow regimes investigated. The cause of 
this would be the secondary flow velocity components of the 
vortices, transporting heated fluid from the concave boundary 
wall inward, while simultaneously displacing the cooler fluid 
particles in the center, causing them to move outwardly 
towards the heated concave wall. 

This study was conducted using a single rectangular cross- 
sectional channel, that incorporated both a straight and a 
curved test section. The results obtained in the straight 
section were compared to those results obtained in the curved 
section, at the same flow velocities, in an effort to 


determine the effects of the Taylor-Gortler vortices on the 
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Peeters OLt feat. Additionally, the results of both. the 
Straight and curved section tests were compared to the 
emeenimental results of Ballard [Ref. 30] and Durao [{Ref. 29}. 
imesceenesults of the straight section that fell into the 
laminar flow region, were also compared to the analytical 
sOlution as given by P. McCuen, et al., [Ref. 35]; while 
those in the transition region were compared to the 
experimental results of Kays and Lueng [Ref. 28], and the 
analytical results of the Dittus-Boelter equation [Ref. 38] 
and Eckert's equation [Ref. 39]. While these equations are 
for fully turbulent flows in circular tubes they have been 
used to compare turbulent flows in other geometry channels 
when the diameter of the tube is replaced by the hydraulic 
diameter (wetted permimeter divided by the cross-sectional 
area). 

The straight section results, for transition flow, were 
also compared to the analytical solutions of Petukhov and 
Pemey |Ret. 25| amd Shibani amd Ozisik [Ref. 26], the former 
for circular tubes while the latter is for parallel plates, 
Do@aewith fully turbulent flows. 

The curved section results for flows in the transition 
region were compared to the experimental results of Brinich 
and Graham [Ref. 31], and also to the solutions of Petukhov 


and Popov, and of Shibani and Ozisik, mentioned above. 


LS 





ieee eR MENTAL WORK 


PeeeweoCRIPTION OF THE APPARATUS 

A channel of rectangular cross-section built as described 
Meeererences 29 and 30, was used to achieve the objectives 
of this study. The channel was manufactured from two 0.635 
Semeameter thick sheets of plexiglas, separated by 0.635 
centimeter spacers whose inside lengths were used as the 
Sides of the channel. 

The channel, as shown in Figure 3, was composed of a 
Straight section, 122.0 centimeters long, followed by a 
curved section of 180 degree arc. The interior concave wall 
of the curved section has a radius of curvature of 30.5 
Seieemmeters, tithe rectangular channel, with an aspect ratio 
of 40, was 0.635 centimeters high by 25.4 centimeters wide, 
with a cross-sectional area of 16.13 Square centimeters and 
a wetted perimeter of 57.07 centimeters. A cross-sectional 
view of the channel, for both straight and curved sections, 
1s shown in Figure 4, 

The working fluid used in this investigation was air, at 
room temperature, which entered the channel via an entrance 
bell constructed of Plexiglas, which was connected to the 
straight section. The entrance bell was designed and 
manufactured in accordance with ASME nozzle standards, with 


an elliptical curved base on a major axis equal to ten inches 
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Miememninoteaxis Of One inch. In an attempt to prevent the 
introduction of dust and other foreign matter into the flow 
Channel, cheese-cloth was attached to the opening of the 
siemencesbell nozzle. After passing through the channel, 
the flow was directed through an aluminum exhaust nozzle 
which was connected, by 3/4 inch PVC piping, to a Fisher and 
POrter Company variable area flow meter model 10A3565A. 

The rotameter, with a Bead Guide Tube number FP-1-35-G-10 
and a 1GSVGT float serial number 68T60, at 100% full scale 
had a flow rate of 0.793 cubic meters of air per minute 
(28.0 standard cubic feet of air per minute at one atmosphere 
pace O° F) . 

The flow of air was drawn through the channel, piping and 
rotameter by a single electrically driven Spencer Turbo 
Compressor, 30 Hp. at 3500 rpm, rated at 550 cubic feet per 
minute at 70°F and one atmosphere. The flow rate was con- 
trolled by the use of valving both at the compressor and at 
Giie exit of the rotameter. 

To obtain the experimental heat transfer data for this 
study, two test sections of the channel were constructed: 

(1) a straight test section, 29.2 centimeters in length, 
located far enough downstream to ensure hydrodynamically 
fully developed flow, with a heated test section area of 
741.9 square centimeters. (2) a curved test section, 28.5 


centimeters in length, subtending an arc of 53.1 degrees, 


ay 





located in the lower half of the concave curved portion of 
Giewelannel, with a test section area of 717.4 square 
Semeimeters. 

In each of the test sections, the outer wall was 
modified as shown in Figures 5 and 6. Temsheet, a carbon 
impregnated porous paper with a property of uniform 
Gmeetrical resistivity, was glued to the interior surface of 
Premouter wall of each Of the two test sections. The flow of 
air was then to be heated by Joulean heating of the Temsheet. 
Pemeetesclectrical resiStamce of Temsheet is not constant, but 
Poammet SLOWLY variable with temperature, a precision resistor 
with an electrical resistance of 2.078 ohms was connected in 
Pemeeco with the Témsheet to allow calculation of the 
instanteous power being supplied. 

The variables to be measured, in this investigation, were: 
@emeeme inlet £Llow temperacure for cach Of the test sections 


CT (A) tne Cutlet Gite weeemperature at the exit of each of 


IN)? 
mimemrest sections (Toya? » (3) the temperature of the heated 
boundary wall of each of the test sections ee: (4) the 
temperature of the unheated boundary wall of each of the test 
sections Ce [Se theweermtcmperaturesat the exit of the 


channel (T (6) the™temperature between the two insula- 


EXIT? ’ 
tion layers for each of the two test sections (Tro) » ee 
Eiemcltave across the precision resistor (Vop) > Cojareue 
heater voltage (Vi), nd Omeicumrow Tate through the 


rOmameter ({Q). 
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iiPrmectpeuaelresmMeasunements required ahove, as well as 
fiimgemer temperature measurements in this investigation, 
were made with Copper-Constantan thermocouples. All the 
thermocouples had been previously calibrated using a 
ROSEMOUNT Communtating Bridge model 920A, and a ROSEMOUNT 
Constant Temperature Bath model 913A, as outlined by Ballard 
[Ref. 30]. 

[weed Of Sixty-four Copper-Constantan, glass insulated, 
30 gauge thermocouples were constructed and inserted in 
e@eet 1c locations to facilitate the measurement of the 
required temperatures. All the thermocouples were wired to a 
Hewlett Packard Data Acquisition System, model 2010c, which 
enabled an automatic printed record of the instantaneous 
reading for each thermocouple location. All the thermocou- 
ples were also connected to a Kaye Instrument Ice Point 
Reference, model K140-4, serial 1780. 

In each of the test sections, five of the thermocouples 
were inserted between the first and second layers of 
ARMAFLEX 22 insulation, and another five inserted between the 
Peeond amd thamd layers. Each of these sets of five 
thermocouples were connected in parallel to facilitate the 
recording of an average value of the temperature (TINS) 
between the respective layers of insulation. Additionally, 
Sem@emot tour thermocouples, also connected in parallel to 


indicate average temperature, were inserted into the channel 


or 








@eethne following three locations; (1) the entrance of the 
meremel, (2) the exit Of the straight test section, and (3) 
the exit of the curved test section. At each of these three 
locations, the thermocouples were spaced five centimeters 
apart in the spanwise direction to allow reading of the 

bulk temperature of the flowing air. 

In each of the two test sections, eight thermocouples 
were placed in direct contact with the Temsheet through 
small diameter holes drilled into the Plexiglas outer walls, 
with the thermocouple beads electrically insulated from the 
Temsheet with ENMAR Heat Resisting Glyceryl Phthalate. The 
final sixteen thermocouples were placed, eight to a test 
section in the Plexiglas unheated boundary wall, opposite 
those positioned in the Temsheet. A detailed sketch of the 
thermocouple placement is provided in Figure 7. The average 
heated outer wall temperature (T ee was provided by those 
thermocouples in contact with the Temsheet. The average 
unheated inner wall temperature (T,;), was provided py tne 
thermocouples embedded in the Plexiglas wall. 

The thermal insulation consisisted of layers ot 
ARMSTRONG ARMAFLEX 22 Sheet Insulation, a flexible foamed 
plastic material, on the heated side of each test section, 
and black foam insulation one-half inches thick elsewhere. 
Three layers of the ARMAFLEX 22 insulation were used, each 


layer had a thickness of 0.636 centimeters, and covered an 
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area slightly larger than the respective test section. 
These layers of insulation were used in the compution of 
heat flux losses to the environment, sample calculations 
of these losses are given in Appendix C. The rest of the 
channel was covered with the one-half inch thick Black 
foam insulation to reduce the losses to the surroundings 
through the plexiglas walls as the heated air transited 
the channel. The insulation was held in place by the use 
of Adhesive Heat Resistant ventilation Duct Tape as seen 
in Figure 8. The entire length of the channel as well as 
all the connections between the channel and the rotameter, 
were sealed with General Electric Silicone Rubber Sealant 
Caulk to ensure that there was no leakage of air into the 
channel or rotameter downstream of the entrance bell 
nezz Lee 

Two. aluminum electrodes, 0.318 centimeters thick, were 
inserted in between the Temsheet and the plexiglas, in each 
GOR Ciicmeesmeeetions, as shown in Figure 4. The plates, 
with the Temsheet and aluminum electrodes were assembled 
i) ttmetommelwaetn plastic terlon screws for the purpose 


of eleéétrical insulation. 


Bo EXG@ERTMENTAL PROCEDURES 
Experiments were conducted with the volumeteric flow 


meeemeanied incrementally from 0.198 to 0.674 cubic meters 
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emai, per Minute to correspond to particular Reynolds 
numbers. The data was recorded and calculations performed 
for each flow rate. The experimental procedures foliowed 
were the same for both the straight and curved test 
sections. The straight test section runs served as a 
Daseline for the curved test section and atlowed tor 
comparison of the results, The appropriate electrical 
power required to heat the Temsheet pvoundary wall to 
approximately 50°U was determined through preliminary 
ewoeraments. This temperature was choosen to allow at 
least a twenty degree celisus temperature difference 
between the unheated and heated boundary wall average 
temperatures and to facilitate achieving steady state 
conditions within a reasonable time span. Estimates of 
time to achieve steady state conditions and decisions as 
to data collection were also based on these preliminary 
test runs. 

The determination as to what constituted a steady 
State condition required careful consideration in view 
of the multiple input variables and their effects on the 
results. The final criteria for steady state was based on 
three input variables: (1) flow rate percentage, as 


recorded from the rotameter; (2) heated boundary wall 
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temperature, as recorded from the thermocouples in that 
area; and (3) the heater voltage, as recorded on the 
mortmeter, When these conditions varied by less that 
two percent over a ten minute period, the experiment 
was considered to have achieved steady state, Under 
these conditions the assumption of steady flow and 
constant heat flux from the heated wall were assumed 
watid. 

Originally, the time required to reach this steady 
State condition was found to be in the range of two 
hours, accordingly all runs were at two and one-half 
hours durations, Data was taken on one hour intervals 
memetie fL£1rSt two hours and then at ten minute intervals 
for the next thirty minutes. Later it was found that the 
higher flow percentages for the curved test section, 
required about thirty minutes more time to full achieve 
memaye state, so test runs with a flow percentage of 
Sixty percent or greater were extended with three 
additional ten minute interval readings for a total time 
Sze ehree hours. 

The instantaneous values for all the temperatures 
were automatically printed by the digital recorder in 
millivolts. A standard thermocouple table was used in 
Scemiuleceron with previously obtained thermocouple 


Etioeition data to convert the millivolts to degrees 


a] 








Setewus. the volumetric flow rate was read directly from 
the installed rotameter as the percent of full scale flow. 
The instantaneous power supplied was calculated using the 


GBelationship: 


Qp = 

RR 
Where the precision resistor voltage (Vpn) and the heater 
voltage (Vis) were recorded in volts from the installed 
voltmeter in the digital Data Acquistion System, and the 
electrical resistance of the precision resistor (Rp, ) had 
Deemepreviously determined to be 2.078 ohms. 

As the area of the curved test section was slightly 
léss than the area of the straight test section, (0.0718 
vs. 0.0742 square meters), the heat flux from the heated 
wall, of the respective test sections was different. The 
room temperature entering the test section was checked 
with the use of a mercury thermometer located in the 
Wlemmpty Of the apparatus. 

During the preliminary test runs of the straight 
section, the thermocouple array recording the exiting 
temperature of the flow of air, was set at three 
millimeters from the heated wall. This was satisfactory 
at the lower volumetric flow rates, but at the higher 


flow rates it led to calculations in which the power 


38 








convected to the air was greater, by approximately fifteen 
to twenty percent, than the power being supplied to the 
heater. It was determined that the thermocouples were 
not indicating the correct bulk temperature of the 
exiting flow of air, but rather a higher value of temp- 
erature, this then resulted in an erroneously high value 
of convected heat transfer, and higher Nusselt numbers. 
To correct this, the system was allowed to run until the 
walls of the downstream curved section achieved steady 
state. Then, calculating the losses, by use of the 
thermocouples in between the insulation layers of the 
curved test section, and ae ae these losses to the 
entire curved section area, the energy loss was calculated. 
The flow exiting the channel was assumed to be thoroughly 
mixed and therefore the temperature as recorded by the 
positioned thermocouples, at the exit of the curved 
section, was used as an accurate representation of the 
average bulk exit temperature. With the final exit 
temperature and the losses now known, the temperature of 
the air as it exited the straight section was calculated. 

Wemcorrect this Condition, the thermocouples were 
inserted further in to the channel, or raised closer to 
the outer wall until the millivolt reading was compatible 
With the calculated correct temperature. After this 


procedure, accomplished at a flow rate of 0.595 cubic 


So 











Samemmecers OL air per minute (7/54 of full scale) was 
completed, additional test runs were conducted at 0.555 
and 0.635 cubic centimeters of air per minute (70% and 
80% respectively), without adjusting the thermocouple 
Peemt1ons, to check the accuracy. In each case, the 
temperature was within 0.25°C of what the energy balance 
indicated it should be. 

It was later discovered that the curved section re- 
Sults had the same problem at flow rates greater then 
Wooo cuDIC Centimeters of air per minute. Since the 
channel construction is such that it does not allow 
further downstream data collection, and since the power 
difference was less than four percent it was decided to 
proceed with the correlation of the data and experimenta- 


tion and note the discrepancies. 
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IV, PRESENTATION OF DATA 


A. ANALYSIS 

imesunitorm electrical resistivity of the Temsheet, 
approximated a constant heat flux surface for the heated 
wall. With insulation on the unheated wall the heat losses 
through that wall were negligible and the wall could be 
considered adiabatic. With the large aspect ratio of 40, 
the channel was considered to be approximately parallel 
plates. The straight portion of the channel upstream of 
the straight test section, was of sufficient length to 
ensure that the flow was hydrodynamically developed for 
the flow velocities investigated. Based on these assump- 
tions, the experimental configuration was modeled as 
forced convection between parallel plates, for hydro- 
dynamically developed and thermally developing flows, 
subjected to a constant heat flux. The boundary conditions 
are one wall at a constant specified heat flux, and the 
other wall insulated. 

Memanalyze this condition, several quanties were 
defined as follows: 

The average heat transfer coefficient between the 


heated wall and the flow of air was defined by the equation: 


Q = hA,, AT 
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where, 'Q,.,.' was the heat convected to the air, ‘h' 
temthe averace heat transfer coefficient, ‘Apr! was the 
area of the heated Temsheet for the section being tested, 
and 'AT' -was the difference between the average heated 
wall temperature ie) and the average bulk temperature 
Seeeene fluid (TR). The average bulk temperature of the 
Fluid was defined as the arithmetic mean of the fluid 
inlet bulk temperature (Tx) and the fluid outlet bulk 
temperature (Toyp) - The delta temperature (AT) could then 


be expressed as: 


- dy 
eee tT = tT = INO 


wo B wo Zz 


The actual heat convected to the air ea was 
calculated separately using th? expression: 
Q =m c el - T ) 
aber Eau OUT IN 
where ‘Spain was the specific heat of the air at constant 


pressure and 'm' was the mass flow rate of the air. The 
mass flow rate, is further defined as the product of the 
volumetric flow rate of the air and the local flow density 
of the air as it transits the rotameter. The volumetric 


MmepweGate Of the air was corrected to standard conditions 


for the rotameter and the local flow density was calculated 
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assuming "perfect gas" behavior at constant pressure and 
Mieeeezing the temperature of the air as it exited the 
channel (Teyr): 


The average Nusselt number was then calculated from the 


average heat transfer coefficient, with the expression: 


coon rc! 
Nu = 
oe 
air 
in which ‘d' is the height ot the channel and eee. 1S hve 


thermal conductivity of the air evaluated at the average 


bulk temperature. 


The Reynolds number was then calculated for each test 


wanes follows: 


Rey = 


Hair ne 


where again 'm! and 'd' are the mass flow rate and height 
of the channel respectively, ae is the cross-sectional 
area of the channel, and ae is the dynamic viscosity 
of the air evaluated at the exit flow temperature. 

For the experiment runs utilizing the curved section, 


the Dean number, defined as: 
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MeeersO cvaitiatced., In this equation, 'd' is again the 
emene ot the channel, and ‘R;' is the radius of curvature 
of the unheated wall surface of the channel (the convex 
wall). 

Mmeeneat losses, by conduction, through the heated 
plate and the insulation, by conduction, were also cal- 
culated for each test section utilizing the expression: 

AT 

NS 
in ° (ex JOS ) A 

INS INS PH 
With AT Fy6' being the difference in temperature, as record- 
ed, between the first and second Layers of AKMAFLEX 22 
insulation, and that recorded between the second and third 
layers. ‘XtNs' was the thickness of each layer of insula- 
einen. ‘Kins! was the thermal conductivity of the 


insulation, and ‘Ap’ was again the area of the heated 


Siaameaac eC . 


In order to verify that Q as calculated by the 


‘air’ 
temperature rise in the air, was actually the convection 
heat transfer from the Temsheet heated surface, an 

estimate of the radiation heat transfer was also calculated 
for each test run. The eight thermocouples in the lower 


Plexiglas surface measured the temperature of that surface 


Baieaitewed the Computation, utilizing the Stefan-Boltzman 
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Waweoft radiation, of the radiated heat transfer (Q..) with 


the expression: 


omeT oo - 
oo WO wi 
iF 
R 


where 'o' is the Stefen-Boltzman constant, and ‘Rp! 1s the 


total resistance to radiation heat transfer between the 


surfaces as computed by: 








In the above equation, ners 1s the emissivity of the 
heated Temsheet wall surface and acer 1s the emissivity of 
the unheated adiabatic wall surface (Plexiglas). 

A sketch of the control volume and the energy balance, 


with sample calculations for one test run, is given in 


Appendix B. 


Eee eoULTS 

The data obtained from the experimental test runs were 
evaluated utilizing the expressions described in the 
previous section, and as illustrated in Appendix B. The 
results of the major variables are shown in tabular form, 
moeiaeiie i for the straight test section, and Table II for 


Boewcurvea test section. A plot of the average Nusselt 
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mummers Verses Reynolds number, is given in Figure 9 for 
Gemparison of the two test sections. The uncertainty bands 
have been indicated for the plotted experimental data and 
eeeemple Calculation for the complete error analysis is 
given in Appendix A. 

As shown, the results indicate an increase in the 
mueemor heat transfer with increasing Reynolds number for 
Memsmmeene Straight and curved test sections. Additionally, 
the heat transfer rate of the curved test section was 
higher than that of the straight test section for each 
Reynolds number investigated. 

In investigating the results obtained, and considering 
those Reynolds numbers below 2000 (1025 based on this 
channel configuration) as laminar, at least squares 


Straight line correlation of the data results in: 


0.67 
Nu = OC es Req 


for laminar flows, and 


— Ooo 
Nu = Oy Rey 


fOn transition regime flows, in the curved test section. 


With correlations of 
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for laminar flows, and 


— Oe Sle 
Nu = 02010 Reg 
peeene transition flows in the straight test section. 

In the laminar flow regime, the increase in heat 
transfer was approximately fifteen percent, while in the 
higher transition flow region the increase was approximately 
thirty percent. This compares favorably to the earlier 
studies by Kreith [Ref. 4], in which he reported that the 
heat transfer rate along concave walls increased from 
twenty-five percent to sixty percent for Reynolds numbers, 


7 and 10°, and 


based on hvdraulic diameter, between 10 
Ballard [Ref. 30], who found an eleven percent increase 


in heat transfer rate for laminar flows. 
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TABLE I 


SU eer eelGHl TEST SECTION RESULTS 





Reg Qair h,. AT Nu 
(W) (W/m“°C) ae) 

784 Ona Zale. Oe, 5 
94] 40.14 Gi 50 ZA O 5 
1099 als ore Lk Zee 5 
Sf 45.13 Zoe 3 25.4 5 
rae es OU 25109 Zo 6 
S71 50.44 26.59 J 57.0 6 
by 9 5a, a0 28.0) Coe 6 
1886 or dL Slee 52 Zod 7 
2042 ao. 24 oa. 99 Zo0 i 
Za0 3 Oul..0:5 Boo 0 24.8 i 
2361 O29 7/ Sis. 5.0 AO 8 
Zo 7 yo) 35.44 Zoe 8 
Zoey 5 69.14 On og 2560 8 
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US tee ek 


Suir ORE CURVED TEST SECTION RESULTS 


Rey De. Gaui h AT 
(W) (W/m2°C) (°C) 

782 114 41.19 oes 24.3 
937 137 AV TUES Posy 24.7 
1092 159 50.88 28.56 24.8 
ne 772 7 53.71 30.49 Aes 
Hea 183 57.16 5 2ew 24.6 
1406 205 58.11 soon Dias 
1570 229 62.63 Siceaye 24.9 
1720 251 69.41 BS 26.0 
1883 as 70.90 38.29 25.8 
2038 297 Ge os 41.66 25.6 
2196 320 80. 36 MOET oe? 
2351 343 83.30 45.98 oe 
2505 366 88.82 46.87 26.4 
2659 388 94.67 49.88 26.4 
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V. DISCUSSION AND CONCLUSIONS 


With the measurement of the unheated wall and the air flow 
temperatures, it was indicated that there was a negligible 
difference between the fluid bulk temperature and the tempera- 
ture of the unheated wall. Additionally, the energy balance 
and sample calculations, provided in Appendix B, demonstrate 
that the radiated heat transfer from the unheated wall sur- 
face is minimal. Thus it was assumed that the heat transfer 
to the flowing air was solely by convection from the heated 
Temsheet. 

The high aspect ratio of the channel as iliscussed pre- 
viously, provided the original basis for the assumption that 
the experimental apparatus, as configured, was one of infinite 
parallel plates. The experimental data, for those Reynolds 
numbers in the laminar flow region, tend to substantiate this 
assumption, in that they approach the theoretical limit for 
average Nusselt number of 5.585, for parallel plates with one 
wall heated at a constant heat flux and the opposite wall 
adiabatic [Ref. 37]. For a Reynolds number, based on hydraulic 
diameter, of 1835, the average value of Nusselt number obtained 
Meas 5.582. 

With the formulation of the problem complete, the 
experimental data obtained can now be compared to the analytical 


solutions and experimental results of the same or similar 


a) dl 





problem. That problem being, flow between infinite parallel 
plates with one wall at constant heat flux and the opposite 
wall adiabatic. 

A comparison of the experimental results obtained in this 
study and the experimental results of Durao [Ref. 29] and 
Ballard [Ref. 30] (who investigated the problem solely for 
laminar flows), are shown in Figure 10, for the straight sec- 
Peuemeeand Figure 11, for the curved section. In both these 
studies, as in this investigation, the channel has an aspect 
ratio of 40, and the experimental procedures were similar. 

It can be seen that the results obtained in this study, in 

the laminar flow regime, appear to be in agreement with the 
data obtained previously for both test sections. Additionally, 
Mumerical correlations of the present data, as given earlier, 
for the laminar flows compares favorably to the numerical 
correlations of both Durao and Ballard, for Reynolds numbers 
greater than 600, as given in Appendix C. In the two investi- 
gations, mentioned above, comparisons were made to the analy- 
E@eeeestudies of McCuen, et. al. [Ref. 35], for heat transfer 
between infinite parallel planes with constant wall temperatures 
and heat flux. Ballard also compared his data to the 
Worse-Schimdt solution [Ref. 36], for heat transfer in tubes 
and annular passages. In both the Durao and Ballard studies, 
the experimental data compared relatively well to the predicted 


PemMeletowOL ) etem, et al. A curve of the predicted solution 
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of the McCuen, et. al. study is shown in Figure 10. The re- 
Strictions of that study were as follows: 

mee the velocity profile is fully established. 

2. the fluid transport properties and density are assumed 
constant. 

3. axial heat conduction is negligible. 

4, viscous energy dissipated is negligible. 

5. conditions are invarient with time (steady state). 

All these restrictions have been accounted for in the 
design of the experimental configuration and/or verified in 
the exnerimental processes. A dimensionless axial length 
coordinate defined as: 


Xx 


Dy Re Pr 


was calculated for 'x' equal to 29.2 centimeters (the length 


Oe 


of the straight test section). In this equation, 'Pr '' is 

the Prandtl number of air, taken as 0.7, ‘Dy! 1s the hydraulic 
diameter (1.239 centimeters), and 'Re ' is the Reynolds number, 
based on the hydraulic diameter. 

Based on the values of x*, it was determined that the 
experiments, for laminar flows, were in the thermal entrance 
region. It was therefore decided to integrate the expression 
of the local Nusselt number as given by McCuen, et. al. [Ref. 


a5]: 


1 


Nu = - 
HSe7 095 xx U/5) - 2(1 +Qwiy xs 


WO 


a0 





for the small values of x, in order to obtain the following 


Seenesston fOr the average Nusselt number: 


2 eel 
enc OME 
wO 


we in (1 


ee) x 
Qwo 


(= (2/3) 
‘ 





In this expression, oe was the heat convected to the fluid 
from the unheated wall, and eS was the heat convected to the 
fluid from the heated wall. In accordance with the previous 
assumptions, Q 4 was considered equal to zero, and or was 
considered equal to Oa? and the average Nusselt number could 
then be calculated, and the expression reduces to: 


a = Oren! Z alee 7o 
Ma = SF In (1 - oerygg ON) 


This value of average Nusselt number was plotted for values 
G@rmesponding to an x* less than 0.042, as this represents the 
maximum value for which the analytic solution for Nusselt 
number is valid. This value of x*, (0.042), for x equal to 
29.2 centimeters, equates to a Reynolds number value of 
approximately 800, based on hydraulic diameter(or 411 based 

On the channel height). 

The graphical representation shows that the experimental 
results of this study as well as those of Ballard and Durao 
plotted above the analytical solution. This difference can 
probably be explained, at least in part, by the difference in 
Peemeertessused in the studies, and by the limitations inherent 
in any experimental work. The side wall effects and the in- 


ability to totally account for all the heat transfer processes 
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amey Or losses could account for the differences of approximately 
me tO cen percent. 

In the comparison of the experimental results obtained in 
this study for those flows not in the laminar regime, the 
comparison and correlation was made to those analytical solu- 
tions and experimental results for completely turbulent flows. 
For the most part, the analytical solutions are valid for 
Reynolds number values greater than 104 (equating to a Reynolds 
memeecr Of 5128 in this study), and tube flow (or other geo- 
metries based on hydraulic diameter). 

A graphical comparison of the results obtained in the 
straight test section runs with the analytical solutions and 
other experimental results is made in Figure 12. The plot 
shows the curve of heat transfer in a straight tube for con- 
stant wall temperature, as given by the Dittus-Boelter equations 
fmet. 38]: 

mo = 0.023 Re oa Pit bak 
and by the empirical relation as given by E.R.G Eckert [Ref. 


os 


_ 0.0396 Re 2°! pr 


OS ae 
Re 
This analytical expressions when extrapolated down show some 
agreement with the values obtained at flow rates between forty 
and seventy percent of the rotameter full scale. Unfortunately, 
those values obtained at the higher flow rates showed less of 


a correlation than those mentioned above. These results, 


Si 








aE a ae 


obtained for volumetric flow rates higher than seventy percent 
may have been affected by the fact that they were near the top 
end of the rotameter scale, and therefore, may inherently be 
somehwat less accurate than the mid-scale readings. The ther- 
mocouple problem discussed previously, may also have contri- 
buted to the decrease in correlation. 

When extrapolated down, the experimental results of Kays 
and Leung [Ref. 28], for heat transfer in annular passages, 
with one wall heated at a constant heat rate, and r* equal to 
1.0, (where r* is the annulus radius ratio r./rv), equivalent 
to parallel plates, appears to correlate well with the present 
results. This extrapolation, was again utilizing a least 
Squares straight line correlation, as given in Appendix C. 
Also shown are the analytical solutions of Petuhkov and Popov 
[Ref. 25], for turbulent heat transfer in round tubes given 


aoe 
wa = (£75) Re. Pr 
k, (£) + k,(Pr)VE7B (Pr 7/9) ~ 1) 
where 'f' the friction factor, has been calculated using the 
Filonenko equation given in Karlekar and Desmond [{Ref. 40] as: 
f = (1.82 log Re - 1.64) 74 
and, 


k, (£) =] + 1.4f 


b mr 
k, (Pr) 7 SPR + i Cay 


iS 
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ieee ttecaecxpressien, given by Shibani and Ozisik [Ref. 


woe teor turbulent heat tmansfer between parallel plates: 


f= a2 "00s 9Re? pr” 
where, 
a = 0.88 pene 
(SoG ts Pr) 
=a ane t 


feo. 53 + 0.5 e 
ween fOr 0.1<Pr <10, and HOMeRe <10°, is also plotted. 

While no conclusions can be drawn from the data obtained 
because it is not in the range of Reynolds numbers of these 
Ewuamcs, it is worthy of note, that the present data as extra- 
polated, and the analytical and experimental works of others 
can be graphically shown to be close. It should also be noted 
Eieetne investigations of Hsu [Ref. 341, Ashton [Ref. 33, and 
Haynes and Ashton [Ref. 32] all plotted experimental results 
between the Petuhkov and Popov expression and that of Shibani 
and Ozisik, as shown graphically in the Haynes and Ashton 
Study above. 

Me@emecurvyed Section data is plotted in Figure 13, and 
compared to the experimental work of Brinich and Graham [Ref. 
31] for flow and mass transfer in a curved channel. The 
analytical solution of Petukhov and Popov, and Shibani and 
Ozisik are also plotted as a reference. A least squares 
[oeneceuwdetoneas given in Appendix C was used to extrapolate 
Mima sea Opesainach and Graham. The accuracy of the data 


points used in this correlation are subject to errors, in that 





the actual values were not given in their study and had to be 
taken from a plot of Stanton number versus Reynolds number as 
Maen im the reference. 

Irrespective of the relatively small differences in the 
meomtts and the correlations, the trends appear to be conclu- 
See ihe heat transfer rates can be seen to increase with 
anereasing Reynolds number, and also that the heat transfer 
Tate in a curved channel is measurably increased over that 
Geen a curved channel 1S measurably increased over that 
rate in a straight channel. This increase in heat transfer 
[mee appears to begin with the development of the Taylor- 
S@ertlcr vortices and continues into the transition and tur- 
bulent flow regions. This result also appears to have been 
reached by K. V. Dement'eva and I. Z. Aronov [Ref. 41], who 
in 1978, investigating hydraulic resistance and heat transfer 
in Curvilinear channels of rectangular cross-section, 
concluded that there was an enhancement of heat transfer 
intensity in curved rectangular corss-section channels as 


compared to straight pipes and channels. 
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VI. RECOMMENDATIONS 


The amount of experimental work remaining in the study of 
heat transfer in curved channels is significant. Additional 
experiments in the transition region and in the fully turbu- 
PemeeeregioOn of fluid flow, should be conducted to refine the 
correlations for average Nusselt number. With this work, it 
would be necessary to design and construct a larger experi- 
mental configuration, capable of higher volumetric flow rate 
and greater temperature differences between the heated surface 
and the fluid inlet temperature. Additionally, care should 
be taken to ensure that the temperature of the exit flow, of 
each test section, is accurately measured. This may be 
accomplished by the use of a mixing cup arrangement at the 
exit flow locations. In a channel designed and constructed 
as the one in this study, care must be taken to allow for the 
removal of the mixing cup at the straight section exit location, 
to ensure that turbulence is not induced in the flow by the 
thermocouples prior to the fluid entering the curved section. 

Also in the way of temperature measurement, a thermocouple 
should be installed at the entrance to the rotameter, or pre- 
ferably some better device or means of determining the 
volumetric flow rate. This would ensure a more accurate 
measurement of the fluid temperature for an evaluation of the 


fluid properties in the mass flow calculation. 
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Cerone) EXPERIMENTAL UNCERTAINTY 


ie Uneertainties for the major variable in the experiments, 
were calculated in accordance by the method described by S. Kline 
and F. McClintoch [Ref. 42]. The estimates of the uncertaintv 
in the measured quantities were made conservatively, and as such, 
Peete 1s Considerable confidence in the uncertainties as calcu- 
lated. An example of the calculations follows. 


The Reynolds number was defined by the equation: 


md 


Reg = i 


- A 
ale Cc 


As such, the uncertainty was calculated by 





The uncertainty in the mass flow rate was determined in a 


Slgjmelar Manner, by its equations: 





yo + YT + 8)? 


dR 
Com 
The uncertainties in the height of the channel, dynamics vis- 

cosity of air, cross-sectional area of the channel, the gas con- 


Seomes pressure, temperature and the reading of the volumetric 


flow rate were obtained from estimates as: 


Gi. dA 
ae). ote, —22= = 00223, —£ = .01976 
d U_ A 
an & Gc 
dR _ aa oe 
S— = .00002, = = .00048, = .00055 
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a= 
0 0200 


The uncertainty in the mass flow rate was then calculated as; 
E~ = 0200 


m 
and the Reynolds number uncertainty was 





dRe, 5 5 > ; 
me = «¥(.0200)* + (.01987)* + (.00223)" + (.01976) “= .0345 
Reg 
Reg = 2042 + 71 


The uncertainty in the Nusselt number was also calculated 


in an identical manner, using the equation: 


= hd 
K 





ar 
and the ee as 





— cl) ae 
ANG 2 yfohy” , dd)? , (airy? 


It was therefore, necessary to first calculate the uncertaintv 


— 


for h. Using the same approach, 





f= Sait _ 
A, aT 
aire 
dh _ ./-dQair,2 PL. 2 dat 2 
Sm ee)? + Gah? + or 
where AT = a - TR 
an 
ms Sei oeeand ss = .0048 
PL 


jie uncertainty in Qair was calculated by: 


Qair = mC (T - T 


Pair Oud IN? 
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and 


Pore given? oo paareze our” Tinje2 

a rs 7 aaron — =a 

Qair m Pair OUR IN 
where 

: lea 

i= 9.0200 , == = 0042, 

m Pair 
and @ our - Tin ) _ oe 

— 2)... ae 2 

out ~TINn 


These calculations result in the following uncertainties: 


Gvair _ 

Sree 0656 
B= .0658 

h 

eee =. 0688 
Nu 


It is therefore, the uncertainty in the difference between 
outlet and inlet temperatures of the flowing air, that is 
ultimately the major source of average Nusselt number uncer- 
tainty. This then is the basis for the recommendation to 
Pewter ascertain the exit temperatures for each test section. 
The calculations, as shown above, were for the 65% straight 
section run. Similar calculations for the curved section at 


65% yield the following; 





dRe, 
Re = ,0345 
d 
Reg = 2038 cae) 
GNU _ 9563 
Nu 


G6 





Na = 9,54 ae US aya ak 
The values for the uncertainties for all the variables in 


the curved section test run are as follows: 


QUANTITY UNCERTAINTY 
Apr .0026 
iN@ .0198 
Cp .0042 
d .0199 
De .0398 
A .0526 
Kair .0003 
Nu .0563 
m .0200 
Q .0200 
Qair BUS 
Re .0345 
Tp BOOS 
Tay melo 
Tout .0106 
ta .0060 
Tour-T in .0482 
AT .0047 
u mOiQe2eZ 
0 .0007 
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Peele Bs oAMPLE CALCULATIONS 


A sketch of the control volume for the energy balance, 
for the straight section, indicating the major heat transfer 
components is shown in Figure 14. The sample calculations 
that follow, show the methods and equations that were used 
EO, Obtain these heat transfer components. The heat transfer 
components were required to satisfy the energy balance, and 
to obtain an average Nusselt number as a function Reynolds 
Mieeoer Of the flowing air. 

iiescurved Section energy balance and calculations are 


Sama ilar 





FIGURE 14. Energy balance in straight section 


68 





SAMPLE CALCULATIONS DATA 


ak 


nak 


Os Ss £D 


PR 


IN 


Twol 
Two2 
Two 3 
Two4 


Twos 
Two6o 
Two 7 
Two 8 


Twil 
Twi2 
Twi3 
Twi4 
Twi5 
Twié 
Twi7 
Twi8 


4 


OUT 


PL 


ROOM 


INS1 
Tins2 


EX 


Dee 


5s 
Oo Senos Vv 


Ser 2 any 
figs 


es ee ee a a NM NM HH YO KY HY HY HY FH OO 


js! 


078 


oO 
~ 524 
OD 
. 840 
840 
oon 
SOs 


eZ 
325 
~245 
2od 


.964 
SOy, 
Ue 
008 


20156 
.047 


04 
ri.00 


.124 


qmele. O76 


mV 
mV 
mV 
mV 
mV 
mV 
mV 


mV 
mV 
mV 
mV 


mV 
mV 
mV 
mV 


mV 
mV 
mV 
mV 


mV 
mV 


.0742 m 


oe 
35 
SS, 


Poe ac 
~54°C 
15> C 


Two = 


Twi = 


Lie ene 
LOR 
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Sumo C 


Bono GC 


80.44°F 








Z 


A. = 0.0016 m 
d = 0.00635 m 
K = 0.02662 w/m°C 
at 
Coair =" 1 0057 - C 
is i: = 29635010 Kg/m.sec 
ay 2 , 
Krys = 1S eee O sn 2 
AX ius = 0.00635 m 
Ee = 0.70 
wo 
ee: = 0.40 
wl 
F ; = 1.0 
WO-W1 =19 2 
0 = 5 eecreen 10 W/m °K 


TEMPERATURE CALCULATIONS 


1. Temperature Difference in Insulation (AT ,,<) 


AT IN = Trys2 ; Trys1 peooato = 55.54 = 6.11°C 


i Sbulk Temperature (T,,) 


3, AN OU. cite. + 27.97 
B 2 2 
3. Mean Temperature Difference (AT) 


AT = T - T, = 50.09 - 24.90 = 25.19°C 
wo B 


POWER CALCULATIONS 


1. Power Supplied (Q_) 
ee eee 


Bee H | (3-.2165)(37.121) 


Qp : ee GOSS = 57.46 W 


2. Heat Lost Through Outer Plate co 


O17 Tas 6.11 


Q = ay 
Lo AXins’*INS**PL —s- (0.00635) /((4.18x10 “) (0.0742)) 


= 2.98 W 
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ss 


4. 


Heat Radiated (Q..) 


a. 


De. 


Heat Convected to Air (Q 


ad. 


D.. 


Radiation Resistance (Rp) 


Rp 


Heat 


Or 








l-e i 
wor. iL es wl 
ApT= wo oy en aa 
= _ += ie) = — 
PL ~wo W1- PL 
Radiated (Qr) 
aaet. «") 
ieee 
Rp 





(5.669x10~°) (0.0742) (323.25 
79 


a, 
4.214 W 


) 


Aaiy: 


4 


SOR) 


Corrected Volumetric Flow Rate (q) 


EG 
il 


r (°R) 
= ex 
QYO = Q ‘fp R 
STANDARD 
L (C0Onmes es) SOR 
60 52969 


Peco Sn eee 


Density (0) 


0 
ii 


RGT=rs, | (@copiecomomenzi7 sem6) 


P 7 101325 
RT 


Wey 7 Ke/m>/sec 
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—_ err 


eee eee ee eeeEeEeEeEEeEee SEs 


c. Mass Flow Rate (m) 


m= qo = (0.0073)(1.177) = 0.0086Kg/sec 
d. Heat Convected to Air (Q.._) 
eed, 


=m Cyaiy Cour - fry? 


iris ) (Le 0057) (2/297-21.32) 


fans 


O20 5519 (KI) sec 


Soro WwW 


AVERAGE HEAT TRANSFER COEFFICIENT (h) 


- _ cair . $3.19 
Ky AT ~ (0.0742) (25.19) 


28.46 w/M?°C 


AVERAGE NUSSELT NUMBER (Nu) 





= hde @smiouconooess) _ 
i —scaiy = T07007) =~ Oe 
REYNOLDS NUMBER (Re) 
_ md _ (0.0086) (0.00635) 
7 ulalA. -5 
aire  (1.983x10 °) (0.0016) 
= 1721.19 
= 1721 


G2 





PENA. CG: CORRELATIONS 


All correlations were obtained using the method of least 
squares (first degree polynomial), as outlined by C. F. Gerald 
[Ref. 43]. In this method the values of Reynolds number (x) 
and Nusselt number (Y) were first converted to the appropriate 
natural logrithemic value, after which the summations of Xs 
= > Ys; and xi ¥5 were calculated. The values of these quanti- 
ties were then placed in a matrix to solve the simultaneous 


equations: 


aN = a4 2X; = 2Y. 
2. 
a 2X; a,2xX) = 2x.Y5 


meme N is the number of data points, and the matrix is 


N 2X5 a, =Y4 
>x 3x2 a ~ Saxe) 
1 i l eel 


After attaining the values of ay and ay» the equation of 


the line became: 


Ln Nu = a_ + a,LnRe 
O A 


OT 
a a 


No =e ° Re 


The quantities for the various studies are given in Table 


LIS 


ae 





The standard deviation for the present study curved section 


morrelation was calculated as: 


2 k =(Y; Z ye) 


al 
. N-n-l 


where: Y. 1s the actual value of average Nusselt number ob- 


tained in the test runs 


ve is the value of average Nusselt number calculated 


from the correlation equation 
N is the number of data points used in the correlation 
n is the degree polynomial used. 

The result of the above method was 


fe = .02426 
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